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Abstract A semi-numerical method is used in order 
to locate the position and calculate the period of peri- 
odic orbits in a 3D composite bisymmetrical potential, 
in a number of resonant cases. The potential consists of 
a 3D harmonic oscillator and a Plummer sphere. The 
outcomes are compared with results found using the 
numerical integration of the equations of motion and 
the agreement is very good. This agreement strongly 
suggests, that semi-numerical methods can be used in 
order to obtain fast and reliable results regarding the 
position and period of the periodic orbits in 3D com- 
posite potentials with a harmonic oscillator part and 
different kinds of perturbations. Comparison with other 
methods of obtaining 3D periodic orbits is discussed. 

Keywords Numerical methods; harmonic oscillators 



1 Introduction 

A large number of studies over the last five decades, 
was devoted to locate the position of periodic orbits in 
a given dynamical system. The basic reason for doing 
this, is that near the stable periodic orbits there are 
similar orbits. On the other hand, orbits starting near 
the unstable periodic orbits can produce large chaotic 
regions. On this basis, there is no doubt, that periodic 
orbits represent the backbone of the whole set of or- 
bits and play an important role in understanding the 
dynamical behavior of a given potential. 

Of particular interest are systems made up of per- 
turbed harmonic oscillators. Systems containing 2D har- 
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monic oscillators were studied by a large number of in- 
vestigators [1-7, 11, 18, 21, 23]. 

In a recent paper [8] , we presented a semi- numerical 
method, which was used in order to find the position 
and the period of periodic orbits in the resonant cases 
uji : u}2 = n : m, in the 2D potential 



^ (x^ + y + a ) 



(1) 



As the results of the above research were in good 
agreement with the outcomes given by the numerical 
integration, we decided to extend our semi-numerical 
method to the 3D potential 



T// \ ^ / 2 2 , 22, 2 2\ 

V[x,y,z} = -[uj^x +uj2y + uj^z ) 



(2) 



where in Eq. (2) wi, uj2 and are the unperturbed 
frequencies of the oscillations along the x, y and z axis 
respectively, while jj. and a are parameters. This poten- 
tial can be considered to describe local motion in the 
central parts of a galaxy. 

The equations of motion for a test particle with a 
unit mass are 
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where the dot indicates derivative with respect to the 
time. 
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The corresponding Hamiltonian is written 



TT ( 2 , 2, 2, 22, 22, 2 2\ 



(a;2 + y2 + ^2 + q;2)1/2 



h, 



(4) 



where Px, Py and Pz are the momenta per unit mass 
conjugate to x, y and z respectively, while ft, > is 
the numerical value of the Hamiltonian, which is con- 
served. The aim of this work is to find the position of 
the periodic orbits and the corresponding period in the 
resonant cases loi : u}2 ■ = n : m : I, where n, m 
and I are positive integers, using a semi- numerical pro- 
cedure and to compare the results with the correspond- 
ing outcomes given by the numerical integration of the 
equations of motion. In our numerical calculations we 
shall use the values ji = 0.001 and a = 0.25, while the 
value of the energy h will be treated as a parameter. For 
the numerical integration of the equations of motion, a 
Bulirsch-Stoer method in double precision is used. The 
accuracy of the calculations is checked by the constancy 
of the energy integral (4), which is conserved up to the 
fifteenth significant decimal figure. 

The paper is organized as follows: In Sciction 2 we 
use some theoretical arguments to find the position 
and period of the 1:1:1 resonant periodic orbits. In the 
same Section, we present a Table for the periods of the 
straight-line periodic orbits. Moreover, we make a com- 
parison between the periods of the straight-line orbits 
derived using the semi-numerical method and the cor- 
responding periods obtained using the numerical inte- 
gration. In Section 3 we compare the semi-numerical 
results with the outcomes given by the numerical inte- 
gration for a number of resonant cases. Finally, a dis- 
cussion and the conclusions of this work are presented 
in Section 4. 



2 The 1:1:1 resonance 

In this Section, we shall study the case when oji = 
W2 = = In this case the system has a spherical 
symmetry and potential (2) takes the form 



V{r) = Ico'r' 



(r2 +a2)^/^^ 



(5) 



y2 _|_ ^2. AH the three components L^, 



where r = x 

Ly and of the test particle's angular momentum and 



the total angular momentum 



Li + Li + Li are 



conserved. The system has a circular periodic orbit. Us- 
ing elementary calculations we find that the radius of 
this circular orbit Tc and the test particle's total angu- 



lar momentum Lc are connected through the following 
equation 

'dv^ 



dr 



while the corresponding value of the energy is 



(6) 



(7) 



As all the three frequencies i^i, W2 and W3 in Eqs. 
(3) are equal the period of the circular orbit is given by 
the formula 

27r 



1/2- 



(8) 



Table 1 A comparison between the period Ta and the period 

r,i found by the numerical integration of the equations of 
motion for the 1:1:1 straight hne periodic orbits when u)\ = 
W2 = = 0.4. 



h 


Ts 


r„ 


0.035 


15.4503 


15.4546 


0.040 


15.4857 


15.4866 


0.04.5 


15.5138 


15.5121 


0.050 


15.5305 


15.5.328 


0.055 


15.5551 


15.5499 


0.060 


15.5706 


15.5644 



Apart from the circular exact periodic orbit the sys- 
tem (5) has an additional type of exact periodic orbit. 
These orbits are straight lines going through the origin. 
We observe that for the motion along the lines 



y = kx, z = Xx, 



(9) 



the equations of motion (3) become identical, while all 

three frequencies of the oscillations are always equal. 
Thus, the straight lines (9) are exact 1:1:1 periodic or- 
bits going through the origin. Our numerical experi- 
ments show that the period of the straight line periodic 
orbits is given by 

27r 



where the frequency Wg is given by the formula 

1/2 



w. 



0.84/x 



(10) 



(11) 



where ra is the radius of the circle V{r) = h. Note that 
all the iso-energetic orbits (9) have the same period. A 
comparison between the period Tg and the period T„ 
found by the numerical integration of the equations of 
motion is made in Table 1 . We see that the agreement is 
very good. In order to quantify the agreement between 
the results obtained from the semi-numerical methods 
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and those derived from the numerical integration, we 
calculate the error T^rr = \{Tn — Ts)/Ts\. For the re- 
sults of the 1:1:1 resonant orbits, presented in Table 1 
we see that Terr < 0.04%. 

The authors would like to clarify that since the na- 
ture and also the shape of the circular and straight-line 
periodic orbits is well-known it is therefore unnecessary 
to provide the corresponding figure 3D plots. 



3 Comparing semi-numerical results to 
numerical outcomes 

In this Section, we shall give semi-numerical formulas 
for the position and the period of the periodic orbits in a 
number of resonant cases. Furthermore, we shall com- 
pare the results given by the semi-numerical method 
with the outcomes derived from the numerical integra- 
tion of the equations of motion. 

We shall start from the resonant case when uti : uj2 ■ 
u)3 = 1:1:2. We look for periodic orbits intersecting the 
z axis at 2: = 2;o, (xq = yo = PzO = 0), with velocities 
Pxo = PyO- Our numerical calculations indicate that the 
value of Pxo is well approximated, if we set 



PxO — 



(12) 



Setting the above initial conditions in the harmonic 
part of the Hamiltonian (4), that is when = 0, and 
solving the resulting equation for zq, we find 

1/2 



2hu}( 



(13) 



1 -I- W3 

Note that Eqs. (12) and (13) give a first approximation 
for the starting point of the 1:1:2 periodic orbit. 

For a better approximation we proceed as follows: 
Using numerical experiments we find that a good ap- 
proximation for the frequency w^s of the 1:1:2 periodic 
orbits is given by the formula 

1/2 

2 , 0.48/i 



W3s = 



{z^ + a^f\ 

Inserting the above value of the frequency wss 
(12) we find 



(14) 



in Eq. 



(15) 



where now Zg is the new value of the starting point 
of the 1:1:2 periodic orbit. Inserting the values of Zg, 

Pxs = Pys, X = y = pz = and = in the 

Hamiltonian (4) we find the following equation 



w 



is 



(zf + a2) 



1/2 



h. 



Solving Eq. (16) we obtain the value of Zg. Remem- 
ber, that the value of Pxs is given by Eq. (15), while the 
value of Pys can be found from the energy integral. The 
period of the 1:1:2 periodic orbits is given by 

47r 

Tg = —. (17) 



Table 2 Starting positions and periods for tlie 1:1:2 periodic 
orbits. Subscript n indicates values found by the numerical 
integration, while subscript s indicates values found using the 
semi-numerical method. The values of the parameters are: 
wi = CJ2 = 0.4 and CJ3 = 0.8. 



h 


Zs 




'pxs 


Pxn 


Ts 


T„ 


0.030 


0.1005 


0.1002 


0.1742 


0.1739 


15.4153 


15.4151 


0.032 


0.1034 


0.1031 


0.1792 


0.1789 


15.4187 


15.4183 


0.034 


0.1062 


0.1059 


0.1842 


0.1840 


15.4221 


15.4218 


0.036 


0.1090 


0.1088 


0.1890 


0.1888 


15.4286 


15.4283 


0.038 


0.1113 


0.1115 


0.1936 


0.1934 


15.4286 


15.4283 


0.040 


0.1142 


0.1141 


0.1982 


0.1981 


15.4318 


15.4317 



Table 2 gives the starting positions and periods of 
the 1:1:2 periodic orbits for several values of the energy 
h. Subscript n indicate values found by the numerical 
integration, while subscript s indicate values found us- 
ing the Eqs. (15), (16) and (17). The values of the pa- 
rameters are: wi = W2 = 0.4 and W3 = 0.8. One can see, 
that the agreement is very good. We calculate the ev- 

VOVSZerr = K^n - Zg) / Zs\, Px,„. = \{Pxn -Pxs)/Pxs\ and 

Terr = \{Tn - Tg)/Tg\. For the results of the 1:1:2 res- 
onant orbits, presented in Table 2 we find that Zerr < 
0.30%, px„, < 0.18%, while Terr < 0.03%. Figure 1 
shows the 1:1:2 periodic orbit, for the above values of 
the parameters when h = 0.034. The initial conditions 
are: zq = z„ = 0.1059, = yo = PzO = 0,PxO = Pxn = 
0.1840, while the value of pyo is found from the energy 
integral in all cases. The integration time for the 1:1:2 
periodic orbit shown in Fig. 1 is t = T„ = 15.4218. 

Now we proceed to study the resonant case when 
wi : W2 : = 2:3:3. We shall study periodic orbits 
intersecting the x axis at a point Xg, (yg = Zg = p^g = 0) 
with velocities p^g = pyg. In this case, the numerical 
calculations suggest that a very good approximation 
for the Xg is given by the formula 



Xg = —J— (/i + O.Ol/x). 

If we insert the a,bove value of with i/g — Zg — Pxs 
and Pzs — Pys in the Hamiltonian (4) we find 



(18) 



= 



Pzs = 

where 

1 



Vg. 



2„2 



(16) Vg = -Ljix 



(x2 + a2)i/2- 



(19) 



(20) 
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Fig. 1 A 1:1:2 resonant periodic orbit starting from the z 
axis. The values of parameters are: = cj2 = 0.4, 0^3 = 
0.8 and h = 0.034. The initial conditions are: zq = Zn = 
0.1059, x-() = yo = p^o = 0, PxO = Pxn = 0.1840, while the 
value of pyo is always found from the energy integral. The 
integration time for this orbits is t = T„ = 15.4218. 



The value of the corresponding period is given by 

^ 47r 
Ts = , 

where 

-1 1/2 



(21) 



Wi + 



,2\3/2 



(22) 



Table 3 Similar to Table 2 but for the 2:3:3 resonant periodic 
orbits. The values of the parameters are: uji = 0)2 = 0.4 and 
U13 = 0.6. 



h 






Pzs 


Pzn 






0.055 


0.7182 


0.7181 


0.1227 


0.1226 


31.0863 


31.0725 


0.060 


0.7501 


0.7500 


0.1275 


0.1274 


31.1221 


31.0970 


0.065 


0.7807 


0.7806 


0.1321 


0.1322 


31.1519 


31.1185 


0.070 


0.8101 


0.8100 


0.1366 


0.1365 


31.1771 


31.1376 


0.075 


0.8386 


0.8385 


0.1410 


0.1411 


31.1984 


31.1546 


0.080 


0.8661 


0.8659 


0.1452 


0.1453 


31.2168 


31.1698 



Table 3 is similar to Table 2 but for the 2:3:3 pe- 
riodic orbits when ui = 0.4 and UI2 = 0J3 = 0.6. As 
one can see, the agreement between the results given 
by the numerical integration and the outcomes from 
the semi-numerical formulas (18), (19) and (17) is very 
good. For the results of the 2:3:3 resonant orbits, given 




Fig. 2 A 2:3:3 resonant periodic orbit, starting from the x 
axis when cji = 0.4, tJ2 = = 0.6 and h = 0.07. The initial 
conditions are: xq = x„ = 0.8100 with yo = 20 = Pxo = and 
PzO = Pzn = 0.1365, while the integration time is t = T„ = 
31.1376. 



in Table 3 we find that Xerr < 0.02%, p,^^^ < 0.07%, 
while Terr < 0.15%. Figure 2 shows the 2:3:3 periodic 
orbit, for the above values of parameters when h = 0.07. 
The initial conditions are: xq = Xn = 0.8100 with 
yo = Zo = PxO = and p^o = Pzn = 0.1365. The in- 
tegration time for the 2:3:3 periodic orbit shown in Fig. 
2 is i = T„ — 31.1376. Here we must emphasize that for- 
mulae (18), (19) and (20) give very good results also for 
the resonant cases loi : LO2 ■ 0J3 = 3:2:2 and Wi : W2 : 1^3 
= 4:3:3. The value of the corresponding period Tg for 
the resonant case 3:2:2 is given by the relation 

T. = (23) 



where 



Wu = 



0.52/Lt 



{xl + a2) 



3/2 



1/2 



(24) 



In the resonant case 4:3:3 the period is given by 
equation 



Stt 

Wu 



where 



Wu 



0.63/x 



{xl + a^f^ 



1/2 



(25) 



(26) 
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Fig. 3 A 3:2:2 resonant periodic orbit, when lui = 0.3, tJ2 = 
LJ3 = 0.2 and h = 0.30. The initial conditions are: xq = x„ = 
1.4908 with yo = zo = PxO = and p^o = Pzn = 0.4478, while 
the integration time is t = Tn = 62.7798. 




Fig. 4 A 4:3:3 resonant periodic orbit, starting from the x 
axis, when uii = 0.8, ui2 = t^'s = 0.6 and h = 0.50. The initial 
conditions are: xo = Xn = 0.7654 with yo = zo = PxO = and 
PzO = Pzn = 0.5600, while the integration time is t = T„ = 
31.3865. 



A comparison between the semi-numerical results 
and results found by the numerical integration is given 



in Tabic 4 for the resonance case 3:2:2 and in Table 5 
for the resonance case 4:3:3. The values of the param- 
eters are: oji = 0.3, a;2 = 0J3 = 0.2 for the Table 4 and 
Wi = 0.8, W2 = = 0.6 in Table 5. In both cases the 
agreement is very good. For the results of the 3:2:2 reso- 
nant orbits given in Table 4, we find that Xerr < 0.03%, 
Pz^„ < 0.02%, while Terr < 0.01%, while for the results 
of the 4:3:3 resonant orbits, presented in Table 5 we see 

that Xerr < 0.01%, p^^^^ < 0.02%, while Terr < 0.003%. 

Figure 3 shows the 3:2:2 periodic orbit, for the above 
values of parameters when h = 0.30. The initial condi- 
tions are: xq = Xn = 1.4908 with yo = zq = p^o = 
and PzO = Pzn = 0.4478. The integration time for the 
orbit shown in Fig. 3 is t = r„ = 62.7798. The 4:3:3 
periodic orbit is given in Figure 4. The values of param- 
eters are as in Table 5 and the value of the energy is 
h = 0.50. The initial conditions are: Xq = Xn = 0.7654 
with yo = zo = PxO = and p^o = Pzn = 0.5600, while 
the integration time is f = T„ = 31.3865. 



Table 4 Similar to Table 2 but for the 3:2:2 resonant periodic 
orbits. The values of the parameters are: u>i = 0.3 and u)2 = 
UJ3 = 0.2. 



h 


X3 




Pzs 


Pzn 


T, 


Tn 


0.15 


1.0541 


1.0544 


0.3177 


0.3176 


62.6896 


62.6952 


0.20 


1.2172 


1.2174 


0.3662 


0.3661 


62.7375 


62.7401 


0.25 


1.3608 


1.3610 


0.4091 


0.4092 


62.7634 


62.7647 


0.30 


1.4907 


1.4908 


0.4479 


0.4478 


62.7794 


62.7798 


0.35 


1.6101 


1.6102 


0.4837 


0.4836 


62.7899 


62.7898 


0.40 


1.7213 


1.7214 


0.5169 


0.5168 


62.7974 


62.7973 



Table 5 Similar to Table 2 but for the 4:3:3 resonant periodic 
orbits. The values of the parameters are: wi = 0.8 and u>2 = 

W3 ^ 0.6, 



h 


Xs 




Pza 


Pzn 


Ts 




0.35 


0.6404 


0.6405 


0.4692 


0.4693 


31.3685 


31.3685 


0.40 


0.6847 


0.6848 


0.5014 


0.5013 


31.3761 


31.3762 


0.45 


0.7262 


0.7263 


0.5315 


0.5316 


31.3819 


31.3820 


0.50 


0.7655 


0.7654 


0.5601 


0.5600 


31.3864 


31.3865 


0.55 


0.8028 


0.8027 


0.5873 


0.5874 


31.3900 


31.3901 


0.60 


0.8385 


0.8386 


0.6133 


0.6132 


31.3929 


31.3930 



Finally, we shall consider three additional resonant 

cases, the cases when Wi : a;2 : = 5:2:2, toi : UI2 ■ 
= 5:3:3 and lui : uj2 '■ <^3 = 5:4:4. We consider periodic 
orbits intersecting the x axis at a point x^, {ys = Zg = 
Pxs = 0) with velocities Pzs = Pys- In all the above 
three cases the numerical calculations suggest that a 
very good approximation for Xs is given by the relation 

Xs = —.l^h. (27) 
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Fig. 5 A 5:2:2 resonant periodic orbit, starting from the x 
ajcis, when u)\ = 0.5, 1^2=1^3= 0.2 and h = 0.40. The initial 
conditions are: xq = Xn = 0.8001 with ya = zq = p^a = and 
P20 = Pzn = 0.5668, while the integration time is t = r„ = 
62.7944. 



Using the above value of Xs with = Zg = Pxs = 0, 
Pzs = Pys in the Hamiltonian (4), we find 



Pzs = \/h-Vs, 
with 

Vs = ^ujfxl - 



(28) 



(29) 



T, = 



(a;2 + a2)i/2" 
The value of the corresponding period is given by 

IOtt 



where 



Wis = 



Wi + 



0.18/i 



{x1 + a^) 



3/2 



1/2 



for the 5:2:2 resonance, while 

1/2 



Wis = 



Wi + 



0.40/x 



(30) 



(31) 



(32) 



(a;2 + a2)3/2_ 

for the 5:3:3 resonance. In the case of the 5:4:4 reso- 
nance the value of the frequency wis is given by the 
formula 

-1 1/2 

(33) 



Wi, 



0.70/i 



(a;2 + a2)3/2^ 

Tables 6, 7 and 8 are similar to Table 5 but for the 
resonant cases 5:2:2, 5:3:3 and 5:4:4 respectively. The 




Fig. 6 Similar to Fig. 5 but for the 5:3:3 resonant case. The 

values of the parameters are: oji = 0.5, uj2 = i^s = 0.3 and h = 
0.40. The initial conditions = 0.9799 with yo = 

20 = Pxo = and p^o = Pzn = 0.5300, while the integration 
time is t = T„ = 62.7837. 




Fig. 7 Similar to Fig. 5 but for the 5:4:4 resonant case. The 
values of the parameters axe: u)i = 0.5, U12 = 0J3 = 0.4 and h = 
0.40. The initial conditions axe: xq = x„ = 1.1313 with yo = 

20 = Pxo = and pzo = pz„ = 0.4907, while the integration 
time is t = T„ = 62.7750. 
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Table 6 Similar to Table 2 but for the 5:2:2 resonant periodic 
orbits. The values of the parameters are: oji = 0.5 and 0J2 = 
W3 = 0.2. 



h 






Pzs 


Pzn 




T„ 


0.35 


0.7483 


0.7484 


0.5303 


0.5304 


62.7858 


62.7867 


0.40 


0.8000 


0.8001 


0..5667 


0.5668 


62.7935 


62.7944 


0.45 


0.8485 


0.8486 


0.6009 


0.6010 


62.7992 


62.8001 


0.50 


0.8944 


0.8945 


0.6333 


0.6334 


62.8036 


62.8045 


0.55 


0.9381 


0.9382 


0.6641 


0.6640 


62.8071 


62.8097 


0.60 


0.9798 


0.9797 


0.6935 


0.6934 


62.8100 


62.8107 



Table 7 Similar to Table 2 but for the 5:3:3 resonant periodic 
orbits. The values of the parameters are: uii = 0.5 and W2 = 
0)3 = 0.3. 



h 






Pzs 


Pzn 






0.35 


0.9165 


0.9166 


0.4960 


0.4961 


62.77.33 


62.7736 


0.40 


0.9798 


0.9799 


0.5301 


().5:«)0 


62.7833 


62.7837 


0.45 


1.0392 


1.0393 


0.5621 


0.5622 


62.7907 


62.7912 


0.50 


1.0954 


1.0955 


0.5923 


0.5924 


62.7964 


62.7969 


0.55 


1.1489 


1.1488 


0.6212 


0.6211 


62.8010 


62.8014 


0.60 


1.2000 


1.2001 


0.6487 


0.6488 


62.8046 


62.8050 



Table 8 Similar to Table 2 but for the 5:4:4 resonant periodic 
orbits. The values of the parameters are: uii = 0.5 and U12 = 
u>3 = 0.4. 



h 






Pzs 


Pzn 


Ts 


T„ 


0.35 


1.0583 


1.0582 


0.4592 


0.4593 


62.7636 


62.7633 


0.40 


1.1314 


1.1313 


0.4908 


0.4907 


62.7754 


62.7750 


0.45 


1.2000 


1.2001 


0.5204 


0.5203 


62.7841 


62.7837 


0.50 


1.2649 


1.2650 


0.5484 


0.5485 


62.7909 


62.7903 


0.55 


1.3266 


1.3267 


0.5751 


0.5750 


62.7961 


62.7956 


0.60 


1.3856 


1.3857 


0.6005 


0.6006 


62.8004 


62.7998 



values of the parameters are: uji = 0.5, W2 = W3 = 0.2 
for Table 6, oji = 0.5, W2 = uj3 = 0.3 for Table 7 and 
oji = 0.5, a;2 = ^3 = 0.4 for Table 8. We see that, in 
all cases, the agreement is very good. For the results of 
the resonant orbits, given in Tables 6, 7 and 8 we find 
that in all < 0.01%, p^^^^ < 0.02%, while 

Terr < 0.001%. Flgiirc 5 shows the 5:2:2 periodic orbit, 
for the above values of the parameters, when h = 0.40. 
The initial conditions are: xq = x„ = 0.8001 with yo = 
Z{) = Pxo = and p^o ~ Pzn = 0.5668. The integration 
time for the 5:2:2 periodic orbit shown in Fig. 5 is t = 
Tn = 62.7944. Figure 6 shows the 5:3:3 periodic orbit 
when h = 0.40, while the values of the parameters are as 
in Table 7. The initial conditions are: xq = x„ = 0.9799 
with yo = 2;o = PxO = and Pzo = Pzn = 0.5300, while 
the integration time is t = r„ = 62.7837. Finally, a 
5:4:4 periodic orbit is shown in Figure 7. The values of 
the parameters are as in in Table 8, while the value of 
the energy is h = 0.40. The initial conditions are: Xo = 
Xn = 1.1313 with yo = zq = p^o = and p^o = p^n = 
0.4907. The integration time for the 5:4:4 periodic orbit 
shown in Fig. 7 is i = T„ = 62.7750. 



4 Discussion and conclusions 

As it was mentioned in Section 1, the periodic orbits 
play an important role in the dynamical behavior of a 
Hamiltonian system. In order to locate periodic orbits 
in 2D and 3D potentials researchers have used a number 
of numerical or analytical methods [9, 10, 12-17, 22, 20]. 

In addition to the above methods there are the semi- 
mimcrical methods. We have presented scmi-nunic!ric;al 
results in a number of previous papers [2-5] for polyno- 
mial potentials. Semi- numerical methods were also suc- 
cessfully used in celestial mechanics [19]. All the above 
strongly suggest, that these methods are a very useful 
tool for the study of the dynamical systems. 

In this paper we have studied the periodic motion 

in a dynamical system composed of a 3D harmonic os- 
cillator and a Plummer potential in several resonant 
cases n : m : I, where n < 5, to < 4 and I < 4. We 
made this choice for two basic reasons: First, because it 
is not possible to study all resonant cases in a given po- 
tential and second, because the important resonances 
are those with values of n, to and I. We considered res- 
onant cases n : m : I in which the periodic orbits are in 
fact two-dimensional, lying on a particular sub-plane of 
the three-dimensional space. This issue is very impor- 
tant and essential for our investigation, since it makes 
easier the construction of the semi-numerical relations 
regarding the position and the period of the periodic 
orbit. However, it is in our future plans, to study and 
provide results about real three-dimensional periodic 
orbits. Here we must note that potential (2) was also 
used by the first author to study local motion in the 
central region of a galaxy [24, 25]. 

The agreement between the results obtained from 
the semi-numerical methods and those derived from 
the numerical integration of the equations of motion 
depends on the value of the energy h. Our numerical 
experiments indicate, that there is a wide range for the 
value of the energy in every resonant case, in which the 
agreement is sufficient enough. However, in every reso- 
nant case we studied in the previous sections, we chose 
six representative values of the energy inside this range, 
in order to present tables with our data. At this point, 
it would be useful to refer to how the numerical results 
indicate or suggest the formulas giving the position and 
the period of the periodic orbits. In other words, to ex- 
plain how the formulas are obtained. This was done on 
a complete empirical basis, using our previous experi- 
ence, simple expressions, which of course, contain the 
main parameters entering the Hamiltonian (4) and give 
outcomes that are close to those given by the numerical 
integration. If the agreement is not satisfactory, we look 
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for a new formula and so on, until we have achieved a 
good result. 

In particular, the aim of this work was to locate 
the position and find the period of the periodic orbits 
in the Hamiltonian (4) in the above mentioned reso- 
nant cases using a semi-numerical procedure and com- 
pare the results with the outcomes derived from numer- 
ical integration. We believe that using semi-numerical 
methods for finding periodic orbits is very important 
in 3D Hamiltonian systems, because in these systems it 
is not easy to find integrals of motion, or to apply an- 
alytical methods. The difficulty increases in the cases 
when we have non polynomial perturbing terms. Fur- 
thermore, the comparison of the results given by the 
semi-numerical method with the results given by nu- 
merical integration shows that our method is fast and 
reliable. On this basis, our semi-numerical method can 
be considered as an effective and sharp tool in order to 
find the position and the period of periodic orbits in 
3D systems. It is in our plans to try to find more com- 
plicated periodic orbits in 3D dynamical systems in the 
near future. 
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